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FOREWORD 


This is a progress report on the research project, "Analysis and Computation of 
Internal Flow Field in a Scramjet Engine," for the period ended January 14, 1991. Special 
attention during this period was directed to investigate "Radiative Interactions in Chemically 
Reacting Supersonic Internal Flows." The work was supported by the NASA Langley 
Research Center (Theoretical Flow Physics Branch of the Fluid Mechanics Division) 
through the grant NAG- 1-423. The grant was monitored by Drs. J.P. Drummond and A. 
Kumar of FLDMD-Theoretical Flow Physics Branch. 
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RADIATIVE INTERACTIONS IN CHEMICALLY REACTING SUPERSONIC 

INTERNAL FLOWS 

S. N. Tiwari and R. Chandrasekhar 
Department of Mechanical Engineering and Mechanics 
Old Dominion University, Norfolk, VA 23529-0247 

ABSTRACT 

The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic 
flows with finite-rate chemistry and radiation, for hydrogen-air systems. The chemistry source 
term in the species equation is treated implicitly to alleviate the stiffness associated with fast 
reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the 
governing equations in time, until convergence is achieved. The specific problem considered is 
the prcmixed flow in a channel with a ten-degree compression ramp. Three different chemistry 
models are used, accounting for increasing number of reactions and participating species. Two 
chemistry models assume nitrogen as inert, while the third model accounts for nitrogen reactions 
and NO x formation. The tangent slab approximation is used in the radiative flux formulation. A 
pseudo-gray model is used to represent the absorption-emission characteristics of the participating 
species. Results obtained for specific conditions indicate that the radiative interactions vary 
substantially, depending on reactions involving HO 2 and NO species, and that this can have 
a significant influence on the flowfield. 
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NOMENCLATURE 


A band absorptance , m -i 
A 0 band width parameter , m -1 
Cj concentration of the j th species , - mole/m 3 

constant pressure specific heat , 

C<» correlation parameter , (N/m 2 ) 'm -1 
£ iota/ internal energy 
e w Planck's function 
fj mass fraction of j th species 

H total enthalpy , J/kg 

h static enthalpy, J/kg 
k thermal conductivity 
Jb> backward rate constant 
kf forward rate constant 
P pressure , N/m 2 
Pj partial pressure of j th species 
qjl total radiative flux 
R gas constant 

S integrated band intensity , (N/m~) m~ “ 

T temperature. K 

u, v velocity in x — and y — directions , m/ 5 
Wj production rate of j th species , kg/m 3 - s 
x, y physical coordinates 
7 ratio of sped f ic heats 
K p Planck mean absorption coefficient 
A second coe fficient of viscosity , wavelength 


vi 



fi dynamic viscosity , kg /m — s 
T) computational coordinates 
p density 

<7 Stefan — Boltzmann constant 
r shear stress 
<j> equivalence ratio 
uj wave number , m -1 



INTRODUCTION 


In recent years there has been a renewed interest in the development of a hypersonic 
transatmospheric aerospace vehicle capable of flying at sub-orbital speeds. A hydrogen-fueled 
supersonic combustion ramjet (scramjet) engine is a strong candidate for propelling such a 
vehicle. For a better understanding of the complex flowfield in different regions of the engine, 
both experimental and computational techniques are employed. Several computer programs have 
been developed [1—4] and applied to gain more insight into the problem involving the flow in 
the various sections of the scramjet module. 

The combustion of hydrogen and air in the scramjet combustor results in absorbing-emitting 
gases such as water vapor and hydroxyl radicals. Existence of such gases makes the study of 
radiation heat transfer an important issue. There are several models available in the literature 
to represent the absorption-emission characteristics of molecular gases [5-10]. One- and two- 
dimensional radiative heat transfer equations for various flow and combustion related problems 
are available [11-19]. In earlier studies [16,18,19], both pseudo-gray and nongray gas models 
were employed to evaluate radiative heat transfer for chemically reacting supersonic flow. Results 
of both models were compared and the pseudo-gray model was found to be computationally 
more efficient. 

Considerable work has been carried out in the past decade to model the chemical kinetic 
mechanism of the hydrogen-air system. A most complete model would involve some 60 
reaction paths [20], rendering numerical solution very difficult, if not impossible. A two- 
step chemistry model, has been used for computing supersonic combustion [4, 16, 18, 19]. 
This model has only four species and two reaction paths, and is useful for preliminary studies. 
However, there are several limitations to this model, such as ignition-phase inaccuracy (i.e. a 
much shorter ignition delay) and also, overprediction of flame temperature and longer reaction 
times. Recent improvements in this area include a 8-species, 14— reaction model [21] and a 
9-species, 18-reaction model [2, 22]. While none of these aforementioned models account 
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for nitrogen reactions (assuming nitrogen as inert), recent developments in this area include a 
15-species, 35-reaction model which accounts for NO x formation and other nitrogen reactions 
in the hydrogen-air system [22]. 

The objectives of the present study are to extend the radiative heat transfer formulation 
used with the global two-step chemistry model [18, 19], to the more complete models namely 
the 9-species, 18-reaction model as well as the 15-species, 35— reaction model. The effect of 
radiative heat transfer in both transverse and streamwise directions is investigated. The finite- 
difference method using the explicit, unsplit MacCormack scheme [23] is used to solve the 
governing equations. 

The flowfield in the combustor is represented by the Navier-Stokes equations and by the 
appropriate species continuity equations [2, 3]. Incorporation of the finite-rate chemistry models 
into the fluid dynamic equations can create a set of stiff differential equations. Stiffness is due to 
a disparity in the time scales of the governing equations. In the time accurate solution, after the 
fast transients have decayed and the solutions are changing slowly, taking a larger time step is 
more efficient. But explicit methods still require small time steps to maintain stability. One way 
around this problem is to use a fully implicit method. However, this requires the inversion of 
a block multi-diagonal system of algebraic equations, which is also computationally expensive. 
The use of a semi-implicit technique, suggested by several investigators [24-26], provides an 
alternative to the above problems. This method treats the source term (which is the cause of the 
stiffness) implicitly, and solves the remaining terms explicitly. 


BASIC GOVERNING EQUATIONS 


The physical model for analyzing the flowfield in a supersonic combustor is described by 
the Navier-Stokes and species continuity equations. For two-dimensional flows, these equations 
are expressed in physical coordinates as. 


dU_ dF dG 

dt dx + dy 


+ H = 0 


0 ) 
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where vectors U, F, G and H are written as. 


P 

pu 

U = I pv 

pE 

L pfj J 

pu 

pu 2 + P + T xx 

F = | puv + T X y 

( pE + p)u + T XX U + T xy V + 7c* + qRz 

pufj - P D^ 
pv 

puv + Tyx 

G = I PV 2 +P+Tyy 

(pE + P)V + T X yV + TyyV + 7 C y + q Ry 

pvfj - pD^ 

0 
0 

H = | 0 

0 

-Wj 


The viscous stress tensors in the F and G terms are given as. 



(2a) 


(2b) 


(2c) 


where A = — \p . The quantities q cx and qcy in the F and G terms are the components of 
the conduction heat flux and are expressed as 


qcx 


qcy 


, dT 

~ k T* 

. dT 

dy 


( 3 ) 
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The molecular viscosity n is evaluated from the Sutherland’s formula. The total internal 
energy E in Eq. (2) is given by 


r> 2 i 2 

E = + 

P 2 


+ S h ]f) 

j=i 


(4) 


Specific relations are needed for the chemistry and radiative flux terms. These are discussed 
in the following sections. 


CHEMISTRY AND THERMODYNAMIC MODELS 
Chemical reaction rate expressions are usually determined by summing the contributions 
from each relevant reaction path to obtain the total rate of change of each species. Each path 
is governed by a law of mass action expression in which the rate constants can be determined 
from a temperature dependent Arrhenius expression. The reaction mechanism is expressed in 
a general form as 


ns ^ ns 

5^7 ijCj f ‘ = l > nr (5) 

)= i j=i 

where ns = number of species and nr = number of reactions. The chemistry source terms wj in 
Eq. (1) are obtained, on a mass basis, by multiplying the molar changes and corresponding 
molecular weight as 


'] = MjCj = Mj (*£■ ~ 7 «'>) 
1=1 

- h. n 


H fi « 


«m | 

m 


> i m 

'm 


m=l 


m=l 

, j = l,na 


( 6 ) 


The reaction rate constants kf , and kb i appearing in Eqs. (5) and (6) are determined from 
an Arrhenius rate expression as 
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where , 


*'• = AiT "' exp {-Wr) 


(7) 


h, = 



( 8 ) 


k 


tq % ~ 



(9) 


The coefficients A , N and E appearing in Eq. (7) are given in Table 1. The term An in 
Eq. (9) denotes the difference in the number of moles of reactants and products. 

The Gibbs energy term in Eq. (9) is calculated as 

ns ns 

AG/?. = *i gi ~ > J = 1 ’ nr (10) 

j=i j=i 


9± 

R 


B, 


= Aj(T - InT) + ^T 2 + 0-T 3 
+ § 7 ’ 4 + § T5 + jPj+Gjr 


( 11 ) 


The gas constant for the mixture is evaluated by a mass-weighted summation over all species 
as 

r = Y.f> R > <12) 

i= i 


P = P RT 


(13) 


RADIATION TRANSFER MODEL 

Evaluation of the energy equation presented in Eq. (1) requires an appropriate expression 
for the radiative flux term, qR . Therefore, a suitable radiative transport model is needed. 
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Various models are available in the literature to represent the absorption-emission characteristics 
of the molecular species [10]. The equations of radiative transport are expressed generally in 
integro-differential forms. The integration involves both the frequency spectrum and physical 
coordinates. In many realistic three-dimensional physical problems, the complexity of the 
radiative transport equations can be reduced by introduction of the tangent-slab approximation. 
This approximation treats the gas layer as a one-dimensional slab in evaluation of the radiative 
flux (Fig. 1). 

Detailed derivations of radiative flux equations for gray as well as nongray radiation have 
been carried out previously [15, 19]. For a multiband gaseous system, the nongray radiative 
flux in the normal direction is expressed as 


<iR{y) = e i - e 2 + 

Ai 


1=1 A 


X 


3u 0 , 


~(y - z) dz + 


(14) 


L 

I 


’ de^{z) 

dz 


Ai 


3 uoi / x 

2T iz - y) 


dz] 


The information on the band absorptance Ai and other quantities is available in the cited 
references. 

For a gray medium, the spectral absorption coefficient k u is independent of the wave 

number, and an expression for the radiative flux is obtained as [5, 16, 19] 

y 

3 f 3w( y — *) 

qR(y) = ej - e 2 + -{ / [e(r) - ei]e" * *dz 
L - ° 

[e(z) — e 2 ]e * k dz) 

y 

It is computationally more efficient to use Eq. (15) in the general energy equation than Eq. (14). 
This is because by differentiating Eq. (15) twice (using the Leibnitz formula) the integrals are 
eliminated and the following inhomogenous ordinary differential equation is obtained : 
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(16) 


1 <P<m{y) 9 _ ,.. s _ 3 My) 

2 j o a $R\y ) i 

k* dy- 4 ac dy 


The solution of Eq. (16) requires two boundary conditions which are given for non-black 
diffuse surfaces as [5] 

'dq R ' 




dy 


= 0 


(17a) 


y=0 


(i-5)[«WU + Tk 


dqR 


= 0 


dy J y=L 

For black surfaces ei — ti — 1 and Eqs. (17) reduce to simpler forms. 


(17b) 


An appropriate model for a gray gas absorption coefficient is required in Eqs. (15) — 
(17). This is represented by the Planck mean absorption coefficient, which is expressed for a 
multi-band system as [5, 19] 

* = */> = <18> 
It should be noted that /cp is a function of the temperature and the partial pressures Pj of 


the species. 


METHOD OF SOLUTION 

The governing equations are transformed from the physical domain (x , y) to a computational 
domain (f , q), using an algebraic grid generation technique similar to the one used by Smith 
and Weigel [27]. In the computational domain, Eq. (1) is expressed as 

(19) 


dU dF dG ~ 

-TT~ + + J— + H = ° 

dt Of OTj 


where 


U = UJ , F = Fy„ - Gx v 
G = Gxt - Fyt , H = HJ 
J = x^y„ - y^Xr, 


( 20 ) 
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Once the temporal discretization has been performed, the resulting system is spatially differenced 
using the explicit, unsplit MacCormack predictor-corrector scheme [23]. This results in a 
spatially and temporally discrete, simultaneous system of equations at each grid point [25, 26]. 
Each simultaneous system is solved, subject to initial and boundary conditions, by using the 
Householder technique [28, 29]. At the supersonic inflow boundary, all flow quantities are 
specified as ffeestream conditions. At the supersonic outflow boundary, non-reflective boundary 
conditions are used, i.e. all flow quantities are extrapolated from interior grid points. The upper 
and lower boundaries are treated as solid walls. This implies a non-slip boundary condition (i.e. 
zero velocities). The wall temperature and pressure are extrapolated from interior grid points. 
Initial conditions are obtained by specifying ffeestream conditions throughout the flowfield. The 
resulting set of equations is marched in time, until convergence is achieved. The details of the 
radiative flux formulation and method of solution are available in [19]. 


RESULTS AND DISCUSSION 

Based on the theory and computational procedures described previously, an algorithm has 
been developed to solve the two-dimensional Navier-Stokes equations for chemically reacting and 
radiating supersonic flows. The extent of radiative heat transfer in supersonic flows undergoing 
hydrogen-air chemical reactions, has been investigated using three chemical kinetics models, 
accounting for increasing number of reactions and participating species. For the temperature 
range considered in this study, the important radiating species are OH, NO and H 2 O. The gray 
gas formulations are based on the Planck mean absorption coefficient which accounts for the 
detailed information on different molecular bands. The radiative fluxes have been computed 
using this ‘pseudo-gray’ formulation. The justification for using this model is provided in [19]. 
Important results for radiative effects of participating species and of chemical interactions were 
presented at the 29 th Aerospace Sciences Meeting in Reno, Nevada, January 7-10, 1991 (AIAA 
Papers 91-0373 and 91-0572). These are provided in this report as Appendices A and B. 
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The specific problem considered is the supersonic flow of premixed hydrogen and air 
(stoichiometric equivalence ratio <f> = 1.0 ) in a channel with a compression comer on the 
lower boundary (Fig. 2). The physical dimensions considered for obtaining results are L = 2 
cm., Xi = 1cm., X 2 = 2cm., Lx = Xi + X 2 = 3 cm. , and a = 10 degrees . The flow is ignited 
by the shock from the compression comer. The inlet conditions which are representative of 
scramjet operating conditions, are Poo =1.0 atm.. Too = 900 K and Moo = 4.0. This same 
flow has been computed by several CFD research groups [4, 18, 19, 21] as a benchmark case. 

Figures 3-7 show the computed results using a 31x31 grid , for temperature and pressure 
as well as H 2 O and OH species mass fractions. Figure 3 shows the flowfield contours for 
temperature and pressure for the non-reacting (cold) case.The oblique shock is seen to arise 
from the compression comer at the lower wall. The hottest regions in the flowfield are in 
the upper and lower boundary layers. Figures 4 and 5 show the effect of the three chemistry 
models on the temperature and pressure profiles, varying along x at the location y = 0.02 
cm from the lower wall (boundary layer region). The temperatures in the boundary layer show 
a gradual increase(Fig. 4). The pressure profiles are plotted at y = 0.13 cm. (inviscid region) 
and show a sharp increase due to the shock (Fig. 5). The ignition-phase inaccuracies of the 
three chemistry models can be seen in Figs. 6 and 7. The shock is occurring after x / Lx = 
0.3. However, the 2-step model predicts ignition before the shock (shorter ignition delay) due 
to the high temperature in the boundary layer. On the other hand, the 18-step model predicts a 
longer ignition delay, at x / Lx = 0.37 (Fig. 6). The 35-step model’s prediction of ignition delay 
appears to be an average of the other two models. Although the three models do not differ much 
in prediction of temperature and pressure profiles, they do differ significantly in predictions of 
species productions (Figs. 6, 7). 

In order to resolve this discrepancy, a grid sensitivity study was carried out to examine 
whether the grid size affects the flow predictions. The results of three grid distributions 1 1 x 
31 , 31 x 31 and 61 x 31 are shown in Fig. 8, and it appears that the 31 x 31 grid is 
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sufficient for the present study. 


The reason for the varying predictions of species production by the three models was further 
examined and the results are shown in Figs. 9 and 10. Figure 9 shows that the Reaction No. 8 
in Table 1 is critical in determining the extent of chemical heat release and H 2 O production. 
Reaction No. 8 deals with production of HO 2 radical. This reaction is absent from the 2-step 
model, while it is common to both the 18-step and 35-step models. Figure 9 shows that the 
35-step model experiences nearly a 30% drop in temperature at the channel exit, when the rate of 
Reaction No. 8 is reduced by a factor of 1000 (effectively cutting off the production of the HO 2 
radical). In contrast, the 18-step model shows a 15% drop in temperature, when subjected to the 
same reduction in rate of Reaction No. 8. This shows that the Reaction No. 8 controls the overall 
H 2 O production occurring in Reaction Nos. 9-18 (Table 1). Due to the high temperatures ( 
" 3000 K) in the flowfield, there is a pool of highly reactive free radicals like H , O, etc. The 
HO 2 radical is converted to the very reactive OH radical, by the free radicals (Reaction Nos. 
11 and 12). This establishes the HO 2 radical as a very important species in promoting flame 
propagation in hydrogen-air flames. A similar study has been carried out in [30]. Since the 
2-step model does not have the HO 2 radical, it predicts lesser amounts of OH and H 2 O. 

It was necessary to determine the reason for the higher sensitivity of the 35-step model to the 
HO 2 radical, as compared to the 18-step model. Figure 10 shows that the Reaction Nos. 21 and 
23 in Table 1 are critical in determining the extent of chemical heat release and H 2 O production. 
Reaction Nos. 21 and 23 deal with production of the NO radical. These reactions are absent 
from the 2-step and 18-step models, whereas they play an important role in the 35-step model. 
Figure 10 shows that the 35-step model undergoes a 30% reduction in temperature, when the 
rates of Reaction Nos. 21 and 23 are reduced by a factor of 1000 (effectively cutting off the 
production of the NO radical). This is nearly the same reduction caused by reducing the rate of 
Reaction No. 8 by a factor of 1000. Due to the high temperatures in the flowfield, the usually 
inert nitrogen dissociates into the highly reactive N free radical. This free radical N is then 
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oxidized in Reaction Nos. 21 and 23, thereby producing the NO radical. This NO radical 
converts the HO 2 radical into the highly reactive OH radical, through Reaction No. 29. This 
confirms that the NO radical is a very important species for flame propagation in a hydrogen- 
fueled supersonic combustor. Since the 35— step model has the NO radical, it predicts higher 
amounts of OH and H 2 O than the 18-step model. 

Figures 11-13 describe flowfield contours affected by exothermic chemical reactions. Com- 
parison of Figs. 3a and 1 1 shows that the boundary layer on the lower wall is considerably 
thickened due to the heat release from chemical reaction. Figures 3b and 12 show that the 
chemical reaction in the boundary layer causes the shock to get curved upwards, thereby in- 
creasing the shock strength. In the case of the 35-step model (Fig. 12b), two weak shocks arise 
from the channel entrance region. These could be responsible for some of the flow instabilities 
observed in the figures. Figure 13 shows the computed H 2 O flame structure for the 2-step and 
35-step models. The 35-step flame(Fig. 13b) appears to be stronger and undergoes a quenching 
effect (discontinuity in the flame sheet) which is not observed in the 2— step flame. 

Based on the above understanding of the chemical kinetics of supersonic hydrogen-air 
flames, the radiative interactions were examined. Figure 14 shows the profiles of the normalized 
stream wise radiative flux qR X predicted by the three chemistry models, along the location y 
= 0.02 cm. from the lower wall. The qRx flux reduces towards the end of the channel due 
to cancellation of fluxes in positive and negative directions. It is seen from Fig. 14 that the 
18-step and 35-step models predict significantly higher amounts (50% more and 100% more , 
respectively) of qR X than the 2-step model. This is because radiative heat tranfer is a strong 
function of temperature, pressure and species concentrations. So the larger values of radiative 
fluxes are caused by higher amounts of H 2 O concentrations, which in turn, depend on reactions 
involving HO 2 and NO species. Figure 15 shows the variations of the normal radiative flux 
qRy along x , at the location y = 0.02 cm. These do not appear to vary significantly between 
the three chemistry models. However, in all three cases, the qRy value increases rapidly after 
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the shock. 


Figures 16-19 show the computed results for reacting flows with and without radiation, for 
the three chemistry models. It is seen from that the 2-step model shows only slight effect of 
radiative interaction, as compared to the 18-step and 35-step models. The 18-step and 35-step 
models predict lower temperature and lower H 2 O and OH concentrations after the shock. 
This is because of the qR X flux, which reduces the total energy. For reacting flows without 
radiation, it was seen from Figs. 6 and 7 that the 18-step model had a longer ignition delay 
(ignition at x / Lx = 0.37), while the 35-step model had a shorter ignition delay (ignition at 
x / Lx = 0.27). Another effect of radiative interactions, seen in Fig. 18, is to nullify this 
difference in predictions of ignition delay. For both 18-step and 35-step models, with radiation, 
the ignition is seen to occur at the same point, x / Lx = 0.33. No such effect is seen on the 
ignition characteristics of the 2-step model. 

Comparison of results in Figs. 16-19 shows that the 35-step model exhibits stronger effect 
of radiative interactions, than the other two models. Accordingly, Figs. 20-22 show flowfield 
contours affected by the 35-step chemistry model, as well as by radiation. It can be seen 
from Figs. 1 lb and 20 that the effect of radiation is to lower the temperature after the shock 
(cooling effect). Figures 12b and 21 show another effect of radiation — elimination of the two 
weak shocks arising from the channel entrance. This is understood from Fig. 14, where the 
qR X radiative flux is higher in the upstream region of the flow. Figures 13b and 22 show that 
the radiative interactions enable the development of a continuous flame sheet, wherein the flame 
quenching effect observed in Fig. 13b is nullified. 

CONCLUSIONS 

The two-dimensional, spatially elliptic Navier-Stokes equations have been used to obtain 
solutions for supersonic flows undergoing finite-rate chemical reactions along with radiative 
interactions. The specific problem considered is of the premixed flow in a channel with a ten- 
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degree compression ramp. The inlet conditions used in the present study correspond to typical 
flow conditions of a scramjet engine. Three different chemistry models were used for parametric 
studies, accounting for increasing number of reactions and participating species. It is seen that the 
radiative interactions vary significantly, depending particularly on chemical reactions involving 
HO 2 and NO species. These reactions have a substantial effect on the flowfield, with regard 
to H 2 O concentration , temperature and pressure. Also, it is observed that the difference in 
the ignition delays of two chemistry models invoving HO 2 reactions is nullified as a result 
of radiative interaction. The results also show that the streamwise radiative flux reduces the 
temperature and concentration of species. This effect is a strong function of the amount of H 2 O 
species concentration. 
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Table 1. Hydrogen- Air Combustion Mechanism [22] 


R PACTION 

Afmoles) 

Nfcm 1 ) 

B(catofle9/gm-molc) 

** following rcnctloru comillnlc the 18- 

step model ** 


(1) O, + H, « OH + OH 

1.70 x 10” 

0 

48150 

(2) Oj ♦ H «-* OH 4 O 

1.42 x 10” 

0 

16400 

0) H, ♦ OH ++ H,() 4 H 

3.16 x 10* 

IS 

3030 

(4) Hi + O «-» OH 4 H 

2.07 x 10” 

0 

13750 

(5) OH 4- OH «-* H 2 0 4 O 

5.50 x 10” 

0 

7000 

(6) H 4 OH 4 M HiO 4 M 

2.21 x 10” 

— 2.0 

0 

(7) H f H [ M h Hj f M 

6.53 x 10” 

—10 

0 

(8) HOi 4 M 

3 20 x III' 1 

— 10 

0 

(9) OH 4 HO, ** Oi 4 H,() 

3.00 x 10” 

0 

1000 

(10) H 4 HOi hH,4 Oj 

2.53 x 10” 

0 

700 

(11) H 4 HOi « OH 4 OH 

1.99 x 10” 

0 

1800 

(12) O 4 HOi 4-* Oi 4 OH 

5.00 x 10” 

0 

1000 

(13) Hf^tHO* 4.0,4 HjO, 

1.99 x 10” 

0 

0 

(14) Hj ♦ HOj 44 H 4 ll]Oi 

3.01 x 10” 

0 

1R700 

(13) OH 4 HjOi 44 HjO 4 HOi 

1.02 x 10” 

0 

1900 

(16) H 4 HiOi 44 HiO 4 OH 

5.00 x 10” 

0 

10000 

(17) O 4 H 2 Oi 44 OH 4 HOi 

1.99 x 10” 

0 

5900 

(IS) lliOi 4 M 44 OH 4 OH 4 M 

1.21 x 10” 

0 

45500 

** remnlning reaction! complete the 35- Hep model •• 


(19) Oi4M«*O404M 

2.75 x 10” 

—1.0 

11 8700 

(20) NirMwNtNiM 

3.70 x 10” 

—1.6 

225000 

(21) N 4 Oi 44 O 4 NO 

6.40 x 10* 

1.0 

6300 

(22) N 4 NO 44 O 4 N} 

1.60 x 10” 

0 

0 

(23) N 4 OH 44 H 4 NO 

6.30 x 10” 

0.5 

0 

(24) H 4 NO 4 M 4* HNO 4 M 

5.40 x 10” 

0 

— 600 

(25) H 4 HNO 44 Hj 4 NO 

4 80 x 10” 

0 

0 

(26) O 4 HNO 44 OH 4 NO 

5.00 x 10” 

0.5 

0 

(27) OH 4 HNO 44 HiO 4 NO 

3.60 x 10” 

0 

0 

(2«) HOi 4 HNO 4. HiO t 4 NO 

2.00 x in” 

0 

0 

(29) HO! 4 NO 44 OH 4 NO* 

3.43 x 10” 

0 

—260 

(30) H 4 NOi 44 OH 4 NO 

3.50 x 10” 

0 

1500 

(3 ») O 4 NO, 44 O* 4 NO 

1.00 X 10” 

0 

600 

(32) NOi 4 M *4 O 4 NO 4 M 

1.16 x 10” 

0 

66000 

(33) M 4 OH 4 NO 44 HNOi 4 M 

5 60 x 10” 

0 

— 1700 

(34) M 4 OH 4 NO* 44 HNO* 4 M 

3.00 x 10” 

0 

—3800 

(35) Oil 4 HNOi 44 HjO 4 NO, 

1.60 x 10” 

0 

0 

** following reactions constitute the global 2 

-step model (4, 16. 

18, 19J *• 

(»"> 111 4 O, 44 2 OH 

11.4 x 10” 

—10.0 

4865 

(2") 2 Oil 4 11* 4* 2 11,0 

2.30 x 10* 4 

—13.0 

42500 
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Geometry for the premixed H2 — Air reacting flow 



Print H 35cold plt || Data from 35-step model 



20 


Fig. 3a Temperature contours for non-reacting (cold) case 
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Temperature profiles (y = 0.02 cm.) 
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Pressure profiles (y * 0.13 cm.) 
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H 2 O mass fraction profiles (y - 0.02 cm.) 
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OH mass fraction profiles (y * 0.02 cm.) 
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Effect of reaction rates 
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Fig. 10 Effect of NO and HO 2 reactions 
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Fig. 11a Temperature contours for reacting (2-step model) case 
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Fig. lib Temperature contours for reacting (35-step model) case 
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Fig. 12a Pressure contours for reacting (2-step model) case 
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Fig. 13a H 2 O mass fraction contours for reacting (2-step model) case 
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Fig. 13b H 2 O mass fraction contours for reacting (35-step model) case 
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Fig. 14 Profiles of stream wise radiative flux (y = 0.02 cm.) 
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Fig. 15 Profiles of normal radiative flux (y ■ 0.02 cm.) 
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Fig. 16b Radiation effects on temperature profiles using 18-step model (y = 0.02 cm.) 
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Fig. 17a Radiation effects on pressure profiles using 2-step model (y 0-13 cm.) 
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Fig. 17b Radiation effects on pressure profiles using 18-step model (y =0.13 cm.) 
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Fig. 17c Radiation effects on pressure profiles using 35-step model (y =0.13 cm.) 
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Fig. 18b Radiation effects on H 2 O mass fractioh profiles using 18-step model (y = 0.02 cm.) 
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Fig. 18c Radiation effects on H 2 O mass 
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Fig. 19a Radiation effects on OH mass fraction profiles using 2-step model (y 0.02 cm.) 
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Fig. 19c Radiation effects on OH mass fraction profiles using 35-step model (y = 0.02 cm.) 
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Abstract 

The two-dimensional, elliptic Navicr-Slokcs equa- 
tions are used to investigate supersonic flows with 
finite-rate chemistry and radiation, for hydrogen-air 
systems. The chemistry source term in the species 
equation is treated implicitly to alleviate the stiffness 
associated with fast reactions. The explicit, unsplit 
MacCormack finite-difference scheme is used to ad- 
vance the governing equations in lime, until conver- 
gence is achieved. The specific problem considered is 
the premixed flow in a channel with a ten-degree com- 
pression ramp. Three different chemistry models are 
used, accounting for increasing number of reactions 
and participating species. Two chemistry models as- 
sume nitrogen as inert, while the third model accounts 
for nitrogen reactions and NO* formation. The tan- 
gent slab approximation is used in the radiative flux 
formulation. A pseudo-gray model is used to repre- 
sent the absorption-emission characteristics of the par- 
ticipating species. Results obtained for specific condi- 
tions indicate that the radiative interactions vary sub- 
stantially, depending on reactions involving H0 2 and 
NO species, and that this can have a significant influ- 
ence on the flowficld. 

Nomenclature 

A band absorptance , m~ 1 
A v band w idlh parameter y m _1 
('j concentration of the j th secies, kg - mole /m 3 
( V conslanl pressure specific heat, J /kg — l\ 

( *o correlation parameter t (N/in J ) l m’ mi 
E total internal energy 
e u , Planck's function 
fj mass fraction of j th species 
H total enthalpy^ J/ktj 
h static enthalpy, J/ky 
k thermal conductivity 
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ki backward rate constant 
kj forward rate constant 
P pressure , N/m 2 
Pj partial pressure of j th species 
qti total radiative flux 
ft gas constant 

S integrated band intensity , (N/m 2 ) ~ 1 m~ - 
T temjterature, l\ 

u, v velocity in x - ami y - directions , m/s 
tij production rate of j th species t kg/m* - s 
x > y physical coordinates 
7 ratio of sj>ecific heats 
K r Planck mean ab.sorpt imt coi f fici* at 
X second coefficient of viscosity , wan h ngth 
p dynamic viscosity , kg/m - s 
£,f/ computational coordinates 
p density 

<x Stefan — Boltzmann constant 
t shear stress 
4> equivalence ratio 
u * wave number , i/i — 1 

Introduction 

In recent years there has been a renewed 
interest in the development of a hypersonic uansat- 
mospheric aerospace vehicle capable of flying at sub- 
orbital speeds. A hydrogen- fueled supersonic combus- 
tion ramjet (scramjci) engine is a strong candidate for 
propelling such a vehicle. For a beucr understanding 
ol the complex flowficld in different regions of the en- 
gine, both experimental and computational techniques 
arc employed. Several computer programs have been 
developed (1-4] and applied to gain more insight into 
the problem involving the flow in the various sections 
of the scramjci module. 

'ilic combustion of hydrogen and air in the scram- 
jci combustor results in absorbing -emitting gases such 
as water vapor and hydroxyl radicals. Existence of 
such gases makes the study of radiation heal transfer an 
important issue. There are several models available in 
the literature to represent the absorption-emission char- 
acteristics of molecular gases (5-10|. One- and two- 
dimensional radiative heat transfer equations for vari- 
ous flow and combustion related problems arc available 
1 1 1- 19|, In earlier studies 1 16,18, 191, both pseudo-gray 
and nongray gas mtxlcls were employed to evaluate ra- 
diative heat transfer for ci»cinically reacting su|KTsonic 
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flow. Results of both models were compared and the 
pseudo-gray model was found to be computationally 
more efficient 

Considerable work has been carried out in the past 
decade to model the chemical kinetic mechanism of the 
hydrogen-air system. A most complete model would 
involve some 60 reaction paths [20], rendering numeri- 
cal solution very difficult, if not impossible. A two-step 
chemistry model, has been used for computing super- 
sonic combustion [4, 16, 18, 19]. This model has only 
four species and two reaction paths, and is useful for 
preliminary studies. However, there are several limita- 
tions to this model, such as ignition-phase inaccuracy 
(i.e. a much shorter ignition delay) and also, overpre- 
diciion of flame temperature and longer reaction times. 
Recent improvements in this area include a 8-specics, 
14-reaction model [21] and a 9-species, 18-reaction 
model [2, 22]. While none of these aforementioned 
models account for nitrogen reactions (assuming ni- 
trogen as inert), recent developments in this area in- 
clude a 15-species, 35-reaction model which accounts 
for NO* formation and other nitrogen reactions in the 
hydrogen-air system [22]. 

The objectives of the present study are to extend 
the radiative heat transfer formulation used with the 
global two-step chemistry model [18, 19], to the more 
complete models namely the 9-species, 18-reaction 
model as well as the 15-species, 35-reaction model. 
The effect of radiative heal transfer in both transverse 
and slreamwise directions is investigated. The finite- 
diffcrcnce method using the explicit, unsplit MacCor- 
mack scheme [23] is used to solve the governing equa- 
tions. 

The flowflcld in the combustor is represented by 
Lite Navier-Stokcs equations and by the appropriate 
species continuity equations [2, 3]. Incorporation of 
die finite-rate chemistry models into the fluid dynamic 
equations can create a set of stiff differential equations. 
Stiffness is due to a disparity in the lime scales of the 
governing equations. In the lime accurate solution, af- 
ter the fast transients have decayed and the solutions 
are changing slowly, taking a larger lime step is more 
efficient. But explicit methods still require small lime 
steps to maintain stability. One way around this prob- 
lem is to use a fully implicit method. However, this 
requires the inversion of a block multi-diagonal system 
of algebraic equations, which is also computationally 
expensive. The use of a semi-implicit technique, sug- 
gested by several investigators [24-26], provides an 
alternative to the above problems. This method treats 
the source term (which is the cause of the stiffness) 
implicitly, and solves the remaining terms explicitly. 

Basic Governing Equations 

The physical model for analyzing the flowflcld in 
a supersonic combustor is described by die Navier- 


Stokes and species continuity equations. For two- 
dimensional flows, these equations are expressed in 
physical coordinates as. 


dU_ dE 

dt + dx + 


90 

ay 


+ n 


o (1) 


where vectors U, F, G and H are written as, 


F = 


G = 


9 

pu 

U — I pv 

pE 

pfji 

pu 

pir + p + r xt 
puv + r r „ 

(pE + p)u + t> x u + r x „v + q cr + <j„ x 
pufj - pD^f 

pv 

puv + r yx 

PIP + P + Tyy 

(pE + p) v + l~XyV ”f* Tyy V “f" ([ cy 4* Q/{y 

p>’f> - pnfy 
r o 
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// = I o 
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L — Wj 


The viscous stress tensors in the F and G terms are 
given as. 



where A = . The quantities q^ and q cy in the 

F and G terms are the components of the conduction 
heat flux and are expressed as 


Hex 


Hey 


-k™ 

dx 

- k d -L 

ay 


(3) 


The molecular viscosity p is evaluated from the 
Sutherland's formula. The total internal energy E in 
Eq. (2) is given by 


P u J + t> 2 

L = + — - — 

P L 


+ V 11 j fj 
j = i 


(4) 
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Specific relations are needed for the chemistry and ra- 
diative flux terms. These are discussed in the following 
sections. 


The gas constant for the mixture is evaluated by 
a mass-weighted summation over all s[>ccics as 


Chemistry and Thermodynamic Model 

Chemical reaction rate expressions are usually de- 
termined by summing the contributions from each rel- 
evant reaction path to obtain the total rate of change 
of each species. Each path is governed by a law of 
mass action expression in which the rate constants can 
be determined from a temperature dependent Arrhenius 
expression. The reaction mechanism is expressed in a 
general form as 

n$ ^ n j 

> 1 = (5) 

;=i j = i 


where ns = number of species and nr = number of 
reactions. The chemistry source terms wj in Eq. (1) 
are obtained, on a mass basis, by multiplying the molar 
changes and corresponding molecular weight as 


w, = 


i = l 

n» 



n$ t 

(7.7 -7.7) 

h . n cl r 


1 


- k n 




, j = 1,118 


( 6 ) 


The reaction rate constants kr • and lq>i appearing 
in Eqs. (S) and (6) are determined from an Arrhenius 
rate expression as 


R = (12) 

P = pRT (13) 


Radiation Transfer Model 

Evaluation of the energy equation presented in Eq. 
(1) requires an appropriate expression for the radiative 
flux term, qR . Therefore, a suitable radiative transport 
model is needed. Various models are available in the 
literature to represent the absorption-emission charac- 
teristics of the molecular species [10]. The equations of 
radiative transport are expressed generally in inicgro- 
diffcrcntial forms. The integration involves both the 
frequency spectrum and physical coordinates. In many 
realistic three-dimensional physical problems, the com- 
plexity of the radiative transport equations can he re- 
duced by introduction of the tangent-slab approxima- 
tion. This approximation treats the gas layer as a one- 
dimensional slab in evaluation of the radiative flux (Fig. 
1 ). 

Detailed derivations of radiative flux equations for 
gray as well as nongray radiation have been carried out 
previously [15, 19], For a multiband gaseous system, 
the nongray radiative flux in the normal direction is 
expressed as 


*/. = AiTN ’ exp {~^) 

wltcre , 


^ ( 8 ) 



Hie coefficients A , N and E appearing in Eq. (7) 
are given in Table 1. The term An in Eq. (9) denotes 
the difference in the number of moles of reactants and 
products. 

The Gibbs energy term in Eq. (9) is calculated as 


— €2 + ^ I / — 
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f r je -.(j)i 4 - r 3 »oi , , 


The information on the band absorpiancc Ai and other 
quantities is available in the cited references. 

For a gray medium, the spectral absorption coef- 
ficient Ku, is independent of the wave number, and 
an expression for the radiative flux is obtained as [5, 
16, 19) 


At.'/,. - - 57 7 ,,iK . j - !.»»• (HI) 

j=i j = l 


Hi 

R 


= Aj(T-lnT)+ ^-T 2 + 0-T 3 
Oi - Ei * 

+ 1 12 ^ + 20 T + Fj + GjT 


( 11 ) 
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It is computationally more efficient to use Eq. (15) 
in the general energy equation than Eq. (14). This is 
because by differentiating Eq. (15) twice (using the 
Leibnitz formula) the integrals are eliminated and the 
following inhomogenous ordinary differential equation 
is obtained : 


i Afl(y) 

/c 2 dy 1 



3 de(y) 
k dy 


06 ) 


The solution of Eq. (16) requires two boundary con- 
ditions which are given for non-black diffuse surfaces 
as [5J 

= o ow 

+ 07b) 

For black surfaces t\ = t 2 ■ 1 and Eqs. (17) reduce 
to simpler forms. 

An appropriate model for a gray gas absorption 
coefficient is required in Eqs. (15) — (17). This is 
represented by the Planck mean absorption coefficient, 
which is expressed for a multi-band system as [5, 19] 

* = 08> 

It should be noted that up is a function of the tem- 
perature and the partial pressures Pj of the species. 


Method of Solution 

The governing equations are transformed from the 
physical domain (x , y) to a computational domain (( , 
»/). using an algebraic grid generation technique similar 
to the one used by Smith and Weigel [27J. In the 
computational domain, Eq. (1) is expressed as 


dU_ dF_ dG 

dt + oz + 


+ H = 0 


where 


0 = UJ , F = Fy n - Gx„ 
G = G Xi - Fy { , // = UJ 
J = X(y„ ~ Vi*,, 


(19) 


( 20 ) 


Once the temporal discretization has been performed, 
the resulting system is spatially differenced using 
the explicit, unsplit MacCormack predictor-corrector 
scheme [23]. This results in a spatially and temporally 
discrete, simultaneous system of equations at each grid 
point [25, 26]. Each simultaneous system is solved, 
subject to initial and boundary conditions, by using 


the Householder technique [28, 29]. At the super- 
sonic inflow boundary, all flow quantities are speci- 
fied as fieestream conditions. At the supersonic out- 
flow boundary, non-rcflcctivc boundary conditions arc 
used, i.e. all flow quantities are extrapolated from in- 
terior grid points. The upper and lower boundaries are 
treated as solid walls. This implies a non-slip boundary 
condition (i.e. zero velocities). The wall temperature 
and pressure are extrapolated from interior grid points. 
Initial conditions are obtained by specifying freestream 
conditions throughout the flowfield. The resulting set 
of equations is marched in time, until convergence is 
achieved. The details of the radiative flux formulation 
and method of solution are available in [19]. 

Results and Discussion 

Based on the theory and computational procedures 
described previously, an algorithm has been developed 
to solve the two-dimensional Navicr-Stokes equations 
for chemically reacting and radiating supersonic flows. 
The extent of radiative heat transfer in supersonic flows 
undergoing hydrogen-air chemical reactions, has been 
investigated using three chemical kinetics models, ac- 
counting for increasing number of reactions and par- 
ticipating species. For the temperature range consid- 
ered in this study, the important radiating species are 
OH and H 2 O. The gray gas formulations arc based 
on the Planck mean absorption cocHicicnt which ac- 
counts for the detailed information on different molec- 
ular bands. The radiative fluxes have been compuicd 
using this "pseudo-gray’ formulation. The justification 
for using this model is provided in [19]. 

The specific problem considered is the supersonic 
flow of premixed hydrogen and air (stoichiometric 
equivalence ratio <f> = 1.0 ) in a channel with a com- 
pression comer on the lower boundary (Fig. 2). The 
physical dimensions considered for obtaining results 
arc L = 2 cm., Xj = 1cm., X 2 = 2cm., L x = X» + 
X 2 * 3 cm. , and u = 10 degrees . The flow is ignited 
by the shock from the compression comer. The inlet 
conditions which arc representative of scramjct oper- 
ating conditions, are * 1.0 atm., T^ = 900 K 
and Moo = 4.0. This same flow has been computed 
by several CFD research groups [4, 18, 19, 21 1 as a 
benchmark case. 

Figures 3-6 show the computed results using a 31 
x 3 1 grid , for temperature and pressure as well as H 2 Q 
and OH species mass fractions, varying along x at the 
location y = 0.02 cm from the lower wall (boundary 
layer region). Figures 3 and 4 show the temperature 
and pressure profiles predicted by the three chemistry 
models. The temperatures in the boundary layer show 
a gradual increase(Fig. 3). The pressure profiles are 
plotted at y = 0.13 cm. (inviscid region) and show a 
sharp increase due to the shock (Fig. 4). The ignition- 
phase inaccuracies of the three chemistry mtxlcls can 
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be seen in Figs. 5 and 6. The shock is occurring after 
x / Lx = 0.3. However, the 2-step model predicts 
ignition before the shock (shorter ignition delay) due 
to the high temperature in the boundary layer. On the 
other hand, the 18-step model predicts a longer ignition 
delay, at x / Lx * 0.37 (Fig. 5). The 35-step model's 
prediction of ignition delay appears to be an average of 
the other two models. Although the three models do not 
differ much in prediction of temperature and pressure 
profiles, they do differ significantly in predictions of 
species productions (Figs. 5, 6). 

In order to resolve this discrepancy, a grid sen- 
sitivity study was carried out to examine whether the 
grid size affects the flow predictions. The results of 
three grid distributions 11 x 31 , 31 x 31 and 61 x 
31 are shown in Fig. 7, and it appears that the 31 x 
31 grid is sufficient for the present study. 

The reason for the varying predictions of species 
production by the three models was further examined 
and the results are shown in Figs. 8 and 9. Figure 8 
shows that the Reaction No. 8 in Table 1 is critical 
in determining the extent of chemical heat release and 
H 2 O production. Reaction No. 8 deals with produc- 
tion of HO 2 radical. This reaction is absent from the 
2-step model, while it is common to both the 18-step 
and 35-step models. Figure 8 shows that the 35-slep 
mode) experiences nearly a 30% drop in temperature 
at the channel exit, when the rate of Reaction No. 8 
is reduced by a factor of 1000 (effectively cuuing off 
the production of the HO 2 radical). In contrast, the 
18-stcp model shows a 15% drop in temperature, when 
subjected to the same reduction in rate of Reaction No. 
8. This shows that the Reaction No. 8 controls the 
overall H 2 O production occurring in Reaction Nos. 
9-18 (Table 1). Due to the high temperatures ( ~ 3000 
K) in the flowfield, there is a pool of highly reactive 
free radicals like H , O, etc. The HO 2 radical is 
converted to the very reactive OH radical, by the free 
radicals (Reaction Nos. 11 and 12). This establishes 
the HO 2 radical as a very important species in pro- 
moting flame propagation in hydrogen-air flames. A 
similar study has been carried out in (30). Since the 
2-step model docs not have the HO 2 radical, it pre- 
dicts lesser amounts of OH and H 2 0. 

It was necessary to determine the reason for the 
higher sensitivity of the 35-step model to the HO 2 rad- 
ical, as compared to the 18-stcp model. Figure 9 shows 
that the Reaction Nos. 21 and 23 in Table 1 are critical 
in determining the extent of chemical heat release and 
H 2 O production. Reaction Nos. 21 and 23 deal with 
production of the NO radical. These reactions are 
absent from the 2-step and 18-step models, whereas 
they play an important role in the 35-slep model. Fig- 
ure 9 shows that the 35-step model undergoes a 30% 
reduction in temperature, when the rates of Reaction 
Nos. 21 and 23 are reduced by a factor of 1000 (ef- 


fectively culling off the production of the NO radical). 
This is nearly the same reduction caused by reducing 
the rate of Reaction No. 8 by a factor of 1000. Due 
to the high temperatures in the flowfield, the usually 
inert nitrogen dissociates into the highly reactive N 
free radical. This free radical N is then oxidized in 
Reaction Nos. 21 and 23, thereby producing the NO 
radical. This NO radical converts the H0 2 radical 
into the highly reactive OH radical, through Reaction 
No. 29. This confirms that the NO radical is a very 
important species for flame propagation in a hydrogen- 
fueled supersonic combustor. Since the 35-slep model 
has the NO radical, it predicts higher amounts of OH 
and H 2 0 than the 18-step model. 

Based on the above understanding of the chemi- 
cal kinetics of supersonic hydrogen-air flames, the ra- 
diative interactions were examined. Figure 10 shows 
the profiles of the normalized streamwise radiative flux 
qiu predicted by the three chemistry models, along 
the location y » 0.02 cm. from the lower wall. The 
q Kx flux reduces towards the end of the channel due 
to cancellation of fluxes in positive and negative direc- 
tions. It is seen from Fig. 10 that the 18-step and 
35-slep models predict significantly higher amounts 
(50% more and 100% more , respectively) of q K * than 
the 2-stcp model. This is because radiative heat tran- 
fer is a strong function of temperature, pressure and 
species concentrations. So the larger values of radiative 
fluxes are caused by higher amounts of H 2 0 concen- 
trations, which in turn, depend on reactions involving 
H0 2 and NO species. 

Figure 1 1 shows the variations of the normal ra- 
diative flux q* y along x , at the location y = 0.02 cm. 
These do not appear to vary significantly between the 
three chemistry models. However, in all three cases, 
the qK, value increases rapidly after the shock. 

Figures 12-15 show the computed results for re- 
ading flows with and without radiation, for the three 
chemistry models. It is seen that the 2-step model 
shows only slight effect of radiative interaction, as 
compared to the 18-step and 35-slep models. The 
18-step and 35-stcp models predict lower tempera- 
ture and lower H 2 0 and OH concentrations after the 
shock. This is because of the q^ flux, which reduces 
the total energy. Comparison of results in Figs. 12-15 
shows that the 35-step model exhibits stronger effect 
of radiative interactions, than the other two models. 

For reacting flows without radiation, it was seen 
from Figs. 5 and 6 that the 1 8-stcp model had a longer 
ignition delay (ignition at x / Lx = 0.37), while the 
35-slep model had a shorter ignition delay (ignition at 
x/Lx* 0.27). Another effect of radiative interactions, 
seen in Fig. 14, is to nullify this difference in predic- 
tions of ignition delay. For both 18-step and 35-step 
models, with radiation, die ignition is seen to occur at 
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the same point, x / L x * 0.33. No such effect is seen 
on the ignition characteristics of the 2-step model. 

Conclusions 

The two-dimensional, spatially elliptic Navier- 
Stokes equations have been used to obtain solutions 
for supersonic flows undergoing finite-rate chemical 
reactions along with radiative interactions. The spe- 
cific problem considered is of the premixed flow in a 
channel with a ten-degree compression ramp. The in- 
let conditions used in the present study correspond to 
typical flow conditions of a scramjet engine. Three 
different chemistry models were used for parametric 
studies, accounting for increasing number of reactions 
and participating species. It is seen that the radiative in- 
teractions vary significantly, depending particularly on 
chemical reactions involving HO 2 and NO species. 
These reactions have a substantial effect on the flow- 
field, with regard to H 2 O concentration , temperature 
and pressure. Also, it is observed that the difference 
in the ignition delays of two chemistry models invov- 
ing HO 2 reactions is nullified as a result of radiative 
interaction. The results also show that the streumwisc 
radiative flux reduces the temperature and concentra- 
tion of species. This effect is a strong function of the 
amount of H 2 O species concentration. 
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Table 1. Hydrogen- A ir Combustion Mechanism (22) 


REACTION 

A (moles) 

N(cm J ) 

B(cak>ries/g m - mole ) 

** following rtaclionj conslilule the 18-slep model *• 

(1) O, + H, «-» OH -f OH 1.70* I0' 5 0 

48150 

(2) O 2 + H <-. OH + O 

1.42 * I0 14 

0 

16400 

(3) H a ♦ OH ♦. HiO . H 

3.16 * I0 7 

1.8 

3030 

(4) Ha + 0«0H + 1I 

2.07 * I0 14 

0 

13750 

(5) OH + OH *« H 2 0 ♦ O 

5.50* 10“ 

0 

7000 

( 6 ) II+OH + Mm HiO + M 

2.21 * 10 11 

—2.0 

0 

(7) II t H t M h H] i M 

6.53 x I0 17 

—1.0 

0 

(8) II + O, + M «-» IIO 2 + M 

3 20 » 10'* 

— 1.0 

0 

(9) Oil ♦ IIOi « Oi + HiO 

5.00 * 10'* 

0 

1000 

(10) 11 ♦ HOi Hi + Oi 

2.53 * 10 11 

0 

700 

(II) H + HOi OH + OH 

1.99 * I0 14 

0 

1800 

(12) O + HOi h Oj * OH 

5 00 * 10'* 

0 

1000 

(13) llOa+IIOi Oi ♦ H 1 O 1 

1 99 * 10“ 

0 

0 

(14) II, -r HO, mH ♦ 11,0, 

301 * 10" 

0 

18700 

(15) Oil + II 1 O 1 w. ll 2 0 ♦ HO} 

1.02 x I0 ' 1 

0 

1900 

(16) II «■ II 1 O 1 «-» H a O ♦ OH 

5.00 * I0 ' 4 

0 

10000 

(17) O ♦ lljOa ~OII ♦ HO] 

1.99 * 10' 7 

0 

5900 

(18) II 1 O 1 *Mh OH OH + M 

1.21 x I0 17 

0 

45500 

** remaining reactions complete the 35 -siep model •• 
(19) 0, + M<-«0 + 0+ M 2 75 x I0 ,f —1.0 

118700 

( 20 ) N, + M«~*N + N + M 

3.70 * 10*’ 

— 1.6 

225000 

(21) N + Oi ~ O * NO 

6.40 * 10’ 

1.0 

6300 

(22) N 4- NO «-* O + N, 

1.60 * I0 ,J 

0 

0 

(23) N ♦ OH *-> H + NO 

6.30 x 10" 

0.5 

0 

(24) II * NO + M UNO ♦ M 

5.40 x 10'* 

0 

—600 

(25) II + UNO Hi + NO 

4 80 x 10'* 

0 

0 

(26) O -f UNO OH ♦ NO 

5.00 x 10" 

05 

0 

(27) OH ♦ UNO w. 11,0 ♦ NO 

3.60 x 10'* 

0 

0 

(28) IIOi + IINO .. 11)0] ♦ NO 

2.00 * I0 11 

0 

0 

(29) IIO, ♦ NO OH + NO, 

3.43 x 10'* 

0 

—260 

(30) II + NO, «-* OH + NO 

3.50 x I0 14 

0 

1500 

(31) O + NO, m O, ♦ NO 

1.00 x 10 '* 

0 

600 

(32) NO, + M *4 O + NO + M 

1.16 x 10 “ 

0 

66000 

(33) M ♦ OH ♦ NO HNOa *■ M 

5 60 x 10 “ 

0 

—1700 

(34) M ♦ OH ♦ NO, ♦-* UNO, ♦ M 

3.00 x 10 “ 

0 

—3800 

(35) Oil ♦ UNO, m 11,0 ♦ NO, 

1.60 x Ml 11 

0 

0 

** following reactions constitute the global 2 -step 

model (4. 16 , 18 . 19 ) •• 

0") II, +0,^2 OH 

11.4 x I0 47 

— 10.0 

4865 

(2 ') 2 OH ♦ H, *4 2 H,0 

2.50 * lO 44 

— 110 

42500 
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Fig. 1 2 Radiation effects on temperature profiles IQ 
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Fig. 13 Radiation effects on pressure profiles (y = 0. 1 3 cm ) 











HO mass fraction H^O mass fraction 















/-^<£ V 

APPENDIX B 

f//4 J /5"3 


AIAA-91-0572 

INVESTIGATION OF RADIATIVE INTERACTIONS 
IN CHEMICALLY REACTING SUPERSONIC 
INTERNAL FLOWS 


S. N. Tiwari, R. Chandrasekhar, and A. M. Thomas 

Old Dominion University 
Norfolk, VA 


J. P. Drummond 

NASA Langley Research Center 
Hampton, VA 


29th Aerospace Sciences Meeting 

January 7-10, 1991/Reno, Nevada 


: or permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 
170 L’Enfant Promenade, S.W., Washington, O.C. 20024 




1 

INVESTIGATION OF RADIATIVE INTERACTIONS IN 
CHEMICALLY REACTING SUPERSONIC INTERNAL FLOWS 

S. N. Tlwari 1 , R. Chandrasekhar 2 , and A.M. Thomas 3 
Old Dominion University, Norfolk, Virginia 23529-0247 

and 

J.P. Drummond 4 

NASA Langley Research Center, Hampton, Virginia 23665-5225 

Abstract 

Analyses and numerical procedures are presented to investigate the radiative interactions 
of absorbing-emitting species in chemically reacting supersonic flow in various ducts. Specific 
attention is directed in investigating the radiative contributions of H 2 O, OH, and NO under 
realistic physical and flow conditions. The radiative interactions in reacting flows are investigated 
by considering the supersonic flow of premixed hydrogen and air in a channel with a compression 
corner at the lower boundary. The results indicate that radiation can have significant influence 
on the flowfield and species production depending on the chemistry model employed. 

Nomenclature 


A band absorptance, m" 1 

A« band width parameter, m’ 1 

Cj concentration of the j lh species, kg-mole/m 3 

toj Planck’s function 

E total internal energy 

fj mass fraction of j 0 * species 

h static enthalpy, J/kg 

k thermal conductivity, J/m-s-k 
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kb 

kf 

P 

Pj 

QR 

R 

S 

T 

u,v 

x,y 

K-UJ 

Kp 

A 

/* 

6 */ 

p 

a 

T 

<t> 


UJ 


backward rate constant 

forward rate constant 

pressure, N/m 2 

partial pressure of j* species 

total radiative heat flux, J/m 2 -s 

spectral radiative heat flux, J/m 3 -s 

gas constant, J/kg-K 

integrated band intensity, (N/m 2 )'^' 2 

temperature, K 

velocity in x- and y-direction, m/s 

production rate of j* species, kg/m 3 -s 

physical coordinates 

spectral absorption coefficient, m 1 

Planck mean absorption coefficient 

second coefficient of viscosity, wave length, m 

dynamic viscosity (laminar flow), kg/m-s 

computational coordinates 

density 

Stefan-Boltzmann constant, erg/S-cm 2 -K 4 
shear stress 
equivalence ratio 
wave number, m 1 


Introduction 

There is a renewed interest in investigating various aspects of radiative energy transfer in 
participating mediums. Radiative interactions become important in many engineering problems 
involving high temperature gases. Recent interest lies in the areas of design of high pressure 
combustion chambers and high enthalpy nozzles, entry and reentry phenomena, hypersonic 
propulsion, and defence oriented research. 
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Basic formulations on radiative energy transfer in participating mediums are available in 
standard references [1-5]. The review articles presented in [6—15] are useful in understanding 
the radiative properties of participating species and the nature of nongray radiation. The validity 
of radiative transfer analyses depends upon the accuracy with which absorption-emission and 
scattering characteristics of participating species are modeled. There are several models available 
to represent the absorption-emission characteristics of molecular species and these are reviewed 
in [12, 13]. These models have been used to investigate radiative interactions in several duct 

flows [16-29]. 

The purpose of this study is to investigate the effect of radiative heat transfer in chemically 
reacting supersonic flow in various ducts under different physical and flow conditions. This pro- 
vides essential information for investigating the effect of radiative interactions in the combustor 
of a supersonic combustion ramjet (scramjet) engine. This hydrogen- fueled engine is proposed 
for propelling transatmospheric hypersonic vehicles. Several basic codes have been developed 
to compute the flowfield in a scramjet engine [21-24]. The combustion of hydrogen and air 
results in absorbing-emitting gases such as H 2 O, OH, and NO. Specific attention, therefore, is 
directed in investigating the radiative contributions of these gases under realistic conditions. In 
essence, the present effort is a continuating of the earlier work conducted in this area of re- 
search [25-27]. Extensive literature survey is provided in the cited references. A comparison of 
different chemistry models used in investigating radiative interactions is presented in [28]. 


Basic Governing Equations 

The physical problem considered to investigate the effect of radiative interactions in super- 
sonic flow are two-dimensional laminar flow between two parallel plates (Fig. la) and within a 
circular tube (Fig. lb). Another geometry is considered to study the effect of shocks and chem- 
ical reactions on the radiative heat transfer and this consists of a channel with a compression- 
expansion ramp (Fig. lc). The governing equations and boundary conditions are provided here 
for all physical problems considered in this study. 
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The physical problem considered for basic understanding of radiative interaction in super- 
sonic flows is two-dimensional variable property laminar flow between two parallel plates. For 
this model, two-dimensional Navier-Stokes equations in fully conservative form are used to 
describe the flow field. These equations, in physical domain, can be written as [27,29] 


M***— 


(1) 


where vectors U, F, G and H are written as 


G = 


P 

pu 

U = | pv 

pE 

ipfj 

pu 

pu 2 + P + T xx 
pUV + T x y 

( pE + p)« + T xx u + T X yV + q cx + mix 

p*fi - 

pv 

PUV + Ty X 
pu 2 +])+ T yy 

(PE + P)V + TyyV + T X yU + q C y + fl/fy 

pvfj ~pD^ 

0 

0 

H = 0 

0 


-Wj 
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The viscous stress terms appearing in the definitions of F and G are given in [27]. The relations 
for conduction heat transfer in x and y directions are given by 


qcx = 




(2) 


The terms q Rx and q Ry represent radiative fluxes in x and y directions, respectively; expressions 
for these are provided in the next section. The total energy flux in a given direction is given by 
the corresponding last term in the definitions of F or G. Consequently, this formulation involves 
all kinds of energy interaction including frictional (aerodynamic) heating. The coefficient 
of viscosity is evaluated by using the Sutherland’s formula and the coefficient of thermal 
conductivity is calculated by using a constant value of the Prandtl number. The total internal 
energy E appearing in U, F, and G is given by 


E = P/p + 


u 2 + v 2 


+ EM 


i= 1 


( 3 ) 


Equation (1) can be used to obtain solutions for all kinds of compressible flows. However, 
boundary conditions and numerical procedures for different flows are quite different. For 
supersonic flows the inflow conditions are specified and outflow conditions are obtained by 
extrapolation. The boundary conditions used along the surfaces are u=o, v=o, dP/dy = o, and 
T=T W . 

The governing equations and boundary conditions for the supersonic flow through a channel 
with a compression-expansion ramp is essentially the same as for the parallel plate geometry. 
However, a strong shock is produced at the compression comer and the flow becomes highly 
reacting from the beginning of the X 2 — coordinate. (Fig. lc). 

The basic governing equations for chemically reacting compressible flow through a circular 

tube can be written as [29,30] 


dU 

dt 


+ 


dV 1 0 , ,,, u/ 

tt + - 7r(y c ) = H/y 

dx y ay 


( 4 ) 



6 


where x and y represent the stfeamwise and radial coordinates, respectively. In Fig. lb, the 
radial coordinate is denoted by r. Thus, both r and y notations are used to represent the radial 
coordinate for the circular geometry. The definitions of vectors U, F, and G in Eq. (4) are same 
as given in Eq. (1) and vector H is expressed as 

H = [0 0 p + Toe 0 y wj] (5) 


For the circular tube geometry, the viscous terms appearing in Eq. (1) are given by 

, / du du u\ du /*_\ 

T — - A (s + * + ih*' , fc <6a) 


' yy 


.( du du v\ dv 

= -\-ai + a y + y)~ '‘dy 


(6b) 


l xy 


( du dv 

= -"U; + aI 


) 


(6c) 




.( du du v\ v 

= —A ( 7 j — (- -x — I — ) ~ 2/x— 
\dx dy yj y 


(6d) 


The boundary conditions for the circular tube geometry are similar to the parallel plate geometry. 


Radiative Transfer Models 


Evaluation of the energy equation presented in Eqs. (1) and (4) requires an appropriate 
expression for the net radiative flux in each direction. A suitable radiative transport model is 
needed to represent the true nature of participating species and transfer processes. In this section, 
a brief discussion of various absorption models is given and essential equations for the radiative 


flux are presented. 
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Absorption Models 

Several models are available in the literature to represent the absorption-emission character- 
istics of molecular species. The total band absorptance of a vibration-rotation band is given by 

A = / [1 - expf/Cu^M^ — w »)l ^ 

Jco 

where is the volumetric absorption coefficient, w is the wave number, is the wave number 
at the band center, X = Py is the pressure path length, and the limits of integration are over the 
entire band pass. Various models are used to obtain the relation for A in Eq. (7). 

The gray gas model is probably the simplest model to employ in radiative transfer analyses. 
In this model, the absorption coefficients is assumed to be independent of frequency, i. e., 
is not a function of w . A convenient model to represent the average absorption coefficient of a 
gray gas is the Planck mean absorption coefficient k v which is defined as [1] 

k p = [°° k* eUT)du/e b (T) (8a) 


For a multiband gaseous system, this is expressed as 

«P = [P J /(-T 4 )]Ee w ,(T)S i (T) (8b) 

where Pj is the partial pressure of jth species in a gaseous mixture, GwiCO is the Planck function 
evaluated at the ith band center, and S*(T) is the integrated band intensity of the ith band. 

As defined in Eq. (8), k p is a property of the medium. When * p is evaluated at the 
temperature of the gas, it is actually a mean emission coefficient and it becomes equal to the 
actual mean absorption coefficient only for the conditions of equilibrium radiation field. For a 
nonuniform temperature field, the mean absorption coefficient used for the optically thin radiation 
is the modified Planck mean absorption coefficient which for black bounding surfaces is defined 

as [1,9] 


f°° 

K m (TT w ) = / k^TMTw)^ 

iJo 


/ efc(Ttp) 


(9a) 
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Note that K m is a function of both the gas temperature and the wall temperature. An approximate 
relation between « p and « m is available for infrared radiation as [ 1 ] 

« m (r,r.) = K,(7 , .)(7'./r) (9b) 

This expression is usually employed in gray gas radiative energy transfer analyses. 

Several other models for the mean absorption coefficient are available in the literature [1,31]. 
Since these models account for detailed spectral information of molecular bands, this approach of 
radiative formulation is referred to as the “pseudo-gray formulation.” The gray gas formulation 
for radiative transport is very useful in parametric studies. 

There are several nongray models available in the literature to represent the absorption- 
emission characteristics of vibration-rotation bands. These are classified generally in four classes, 
(1) line-by-line (LBL) models, (2) narrow band models, (3) wide band models, and (4) band 
model correlations. A complete discussion on usefulness and application of these models is 
provided in [12, 13]. For many engineering applications, wide band model correlations provide 
quite accurate results. The most commonly used wide band model correlations are due to 
Edwards [5, 11) and Tien and Lowder [9]. The Tien and Lowder correlation for the total band 
absorptance is a continuous correlation and is given by the relation 

A(u,/?)= A(u, /3)/A 0 = ln{uf(t)[(u + 2)/(u + 2f(t))] + 1} (10) 

where 

f(t) = 2.94(1 - exp(— 2.60t)],t = >3/2 

and u = SX/Ao is the nondimensional path length, fi = 2irj/d is the line structure parameter, 7 
is the line half width, S is the integrated band intensity, and A 0 is the band width parameter. 
Equation (10) provides accurate results for pressures higher than 0.5 atmosphere [12, 13]. 

Spectral properties and correlation quantities for various radiation participating species are 
available in [5, 9, 11]. These are useful in gray as well as nongray radiative formulations. 
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Radiative Flux Equations 

For many engineering and astrophysical applications, the radiative transfer equations are 
formulated assuming one-dimensional planar systems. For diffuse nonretlecting boundaries and 
in absence of scattering, the expression for the spectral radiative flux in the normal direction 
is given by [1, 19] 


qR w (y) = eiw - e 2 c 


yy 

' 3 ‘ 

/ Fi w (z)« a ,exp 

z y ) 

/ O 


f L 

' 3 , 

/ F 2 u ,(z)/c w exp 

-^-y) 

J y 



( 11 ) 


where 


F w (z) = e w (z) - ei w> F 2a ,(z) - e 2 (z) - e 2u > 

It should be pointed out that the exponential kernel approximation has been used in obtaining 
Eq. (11). The total radiative flux in a given direction is expressed as 

( 12 ) 


<IR 


-f 


<1 llu, 


A combination of Eqs. (11) and (12) provides a proper form of total radiative flux equation 
for obtaining nongray solutions of molecular species. Any convenient absorption model can be 
used to obtain nongray results. 

For a gray medium, Eq. (11) reduces to a simpler form and upon differentiating the resulting 
equation twice, the integrals are eliminated and there is obtained a nonhomogeneous ordinary 
differential equation as |1, 16, 27 1 

1 d 2 qtt(y) _ — qj^(y) = — — — — — (13) 


dy 2 


4 


K dy 


where k=k p . Equation (13) is a second order differential equation and, therefore, requires two 
boundary conditions. For nonblack diffuse surfaces, these are given as 

fd, 


(^-i)ha (y )U-2: 


1 <|R 

L dy 


= 0 


y=° 


(14a) 
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J qR 


dy 


= 0 


(14b) 


y=L 


Equation (13) along with boundary conditions can be used to obtain the energy equation for 
gray gas radiative interaction. For black walls and Ti = T2, the boundary conditions for Eq. 
(13) become 

qit(l/2) = 0, ^q R (0) = — (dq R /d£)( =0 > T 0 = nL,(=y/L (15) 

z t q 


For a black circular tube, the spectral radiative heat flux in the radial direction is given by 
the expression [17] 

a r*l 2 r r r r 1 

dr ; 

cos 7 ' ' ] 

b 


qR - (r) {£„, 


r° -1 / 

— I F u ,(r / ) «u,aexp 

+ / F w (r') k w aex 

J r sin 7 


[ cos 7 

b/Cy 


(r' - r) 


dr 


(16) 


cos 7 


(r -F r* — 2r sin 7)dr ; 


d7 


where F u ,(r')=e u ,(r') - eu/(T w ) and constants a and b have values of unity and 5/4, respectively. 
A combination of Eqs. (12) and (16) provides a convenient form of the total radiative flux for 
nongray analyses. 

For a gray medium, the expression for the total radiative flux can be obtained from differential 
approximation as [1, 17] 


d_ 

dr 


1 d 


9 dT 4 

- T^qil = — 

4 dr 


For a black tube, the boundary conditions for Eq. (16) are found to be 

, qn(0) = 0 , r G = ki 0 ,£ = r/r 0 


?qit(i) = “ 


7:37(^11) 


bUd£ 


(17) 


(18) 


f=i 


Equation (17) is used along with Eq. (18) for general one-dimensional gray gas formulation 
and analyses. 

With certain modifications, the radiative flux equations presented in this section can also 
be used to investigate the radiative interactions in the flow direction. The procedure for doing 
this is provided in [27]. 
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Chemistry and Thermodynamic Models 


Chemical reaction rate expressions are usually determined by summing the contributions 


from each relevant reaction route (or path) to obtain the total rate of change of each species. 
Each path is governed by a law of mass action expression in which the rate constants can be 
determined from a temperature dependent Arrehenius expression. The reaction mechanism is 


expressed in a general form as 



( 19 ) 


where ns= number of species and nr= number of reactions. The chemistry source terms in Eqs. 
(1) and (4) are obtained, on a mass basis, by multiplying the molar changes and corresponding 
molecular weight as 






h n c m m - n 


( 20 ) 


m= 1 


m—\ 


The reaction constants and appearing in Eqs. (19) and (20) are determined from an 
Arrhenius rate equation given by 

k ti = Ai T Ni cxp(-|L) (21) 


where 

( 1 \ AN / — AGr;\ 

^bi = ^fi/^eqi t ^eqi — ^ J ex P ^ jpp J 

The coefficients A, N, and E appearing in Eq. (21) are given in Table 1 and the Gibbs energy 
term A Gm is calculated as 


ns ns 

AG Ri = Tij© - W ’ j = nr 
j=l j=l 


( 22 ) 


| = Aj(T - In T) + ^ T 2 + ^ T 3 
+^ T ‘ + | T5 + F j + Gj T 


where 
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The gas constant for the mixture is evaluated by a mass-weighted summation over all species 
as 

R = £ fj R j (23) 

j=i 

The equation of state for the mixture is written as 


P = p RT 


(24) 


Method of Solution 


The governing equations are transformed from the physical domain (x, y) to a computational 
domain (£, rj) using an algebraic grid generation technique similar to the one used by Smith and 
Weigel [32], The grid spacing is kept uniform in the flow direction and compressed near the 
boundaries in the normal direction. The governing equations, Eqs. (1) and (4), are expressed 
respectively in the computational domain as 



dP 

d( 


OG - T n 
+ — + H = 0 

or) 


(25) 


where 



dP ]_d_ 

+ V 


(vG) = H/r? 


U = UJ, F = Fy v -Gx,,G = Gx^ - Fy f , 


(26) 


H = HJ , J = y v - x. 

The temporal descretization procedure for Eqs. (25) and (26) in given in [25-27]. Once 
this has been performed, the resulting system is spatially differenced using the explicit unsplit 
MacCormack predictor-corrector scheme [33]. This results in a spatially and temporally discrete, 
simultaneous system of equations at each grid point [34,35]. Each simultaneous system is solved 
using the Householder technique [36,37], and is marched in time until convergence is achieved. 
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The details of gray as well as nongray radiative flux formulations and solution procedure 
available in [19,27]. For the gray gas model, the governing equations are discretized using 
a central difference scheme. The descretization of Eqs. (13) and (17) results, respectively, in 


are 


qj-i 




LBjVCj Ay,/ -I 


+ 7«j 


9j + 


Aj/*j 


Mj+i 


= mis 


(28) 


2 

7 — H 


Aj 2yjCj 


qj+i 


KU 


Bj VCj A yj 


9 o 

+ 4 K j 


where 


J_ _J 1 

+ y I + 2yjCj 2yjAyj 


qi + 


Aj/?j ‘iyjAvj 


qj-l = aHS 


Aj = Ayf(l + ft\) , 0j = Ayj(l + ft j) , C'j — /^j Ayj , 


RIIS = 1.5/ci 


e j+i ~ e j , e j ~ e i-i 
[ Ayj Ayj 


, Ayj = yj - yj-i. PJ- — 


yj+i - yj 


yj - yj-i 


(29) 


Equations (27) and (28) along with boundary conditions given by Eqs. (15) and (18) form 
tridiagonal systems of equations which can be solved efficiently by the Thomas algorithm. 

In the nongray gas formulation, the divergence of the radiative fiux is evaluated using a 
central differencing scheme and is treated as radiative source term in the energy equation. Since 
the radiative flux terms involves integral formulation, unlike other flux terms which are only u] 
a differential form, it is uncoupled and treated separately. 


Physical Condit ions and Data Source 

The physical conditions for which specific llowfield analyses and computations are needed 
are discussed in [2 1-28 1. In this work selected parametric studies have been conducted for 
certain flow and physical conditions. Radiation participating species considered are H 2 0, OH, 
and NO. Radiative properties of these species are available in [5, 9, 11-13]. Different amounts 
of these gases, in combination with air, are considered for parametric studies. Essential data for 
the chemistry model employed are obtained from Refs. 38-40 and these are provided in Table 1. 
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For basic studies, the physical dimensions considered for the channel are L=3 cm and L x — 10 
cm, and for the circular tube they are L=D=3 cm and L»= 10 cm. Certain results, however, 
were obtained also for other dimensions. 

Results and Discussion 

Based on the theory and computational procedure described in the previous sections, an 
existing computer code was modified to solve the two-dimensional Navier-Stokes equations 
for radiating supersonic laminar flows between two parallel black plates. A similar code was 
developed for radiating supersonic flows in a circular tube. In most cases, the radiative interaction 
was considered only in the normal direction. Extensive results have been obtained for pure H 2 0, 
OH, and NO as homogeneous participating species, and for different mixtures of these species 
with air. Selected results are presented and discussed in this section. * _ 

For the parallel plate geometry (3 cm x 10 cm), a comparison of the divergence of radiative 
flux for general (nongray), gray, and their optically thin limit models is presented in Fig. 2 for 
two different y-locations (y=0.2 and 1.5 cm). The inflow conditions for this case are Poo = 1 atm, 
= 1.700 K, Moo = 4.3, f Hj o = 0.5, fo 2 = 0.1, and f N , = 0.4. The gray formulation is based 
on the modified Planck mean absorption coefficient which accounts for the detailed information 
on different molecular bands. The range of optical thickness calculated in [27] was found to 
be between 0.0003 and 0.4. Thus, for the physical model and inflow conditions considered, 
the radiative interaction is essentially in the optically thin range. No significant difference in 
results is observed for the two y-locations. This is a typical characteristic of the optically thin 
radiation [1,10]. The solution of the gray formulation requires about ten times less computational 
resources in comparison to the solution of the nongray formulation [27]. 

As mentioned earlier, the Planck mean absorption coefficient (or /c m ) is considered to 
be an optically thin radiation absorption coefficient, although it has been used in other optical 
ranges as well [1,9]. The appropriate absorption coefficient for the optically thick radiation is the 
Rosseland mean absorption coefficient kr. It has been pointed out in [1] that if the medium is 
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gray then kr=kp=k; otherwise «p>/CR. Thus, use of /cp (or Km) in pseudo-gray gas formulation 
will provide maximum influence of radiative interaction. The k p values for H 2 O, OH, and NO 
have been calculated from Eq. (8b) by employing radiative properties available in [5, 9, 11] 
and these are illustrated in Fig. (3). Figure 3a shows the results of 1(X)% homogeneous species 
whereas results in Fig. 3b are for different mixtures. The results provide indication of radiative 
ability of different species at a given temperature. Values of k p for other species are available 
in [1, 9]. 

The results of supersonic entrance region flow between parallel plates are presented in Figs. 
4-10 for different physical and inflow conditions. As mentioned earlier, the basic physical 
dimensions considered for the channel are L=3 cm and L x =10 cm. In most cases, the inflow 
conditions considered are Poo = l atm, Too=l,700 K, Uoo = 2574 m/s (Moo = 3.0) with varying 
amounts of radiation participating species in combination with air. Certain variations in physical 
and inflow conditions are also considered for parametric studies. The chemical reactions are 
not considered in obtaining the results of Figs. 4-10, and radiative flux results are presented 
only for the normal direction. 

Results for radiative flux are illustrated in Figs. 4 and 5 as a function of the nondimensional 
y-coordinate. For P=1 atm, the results presented in Fig. 4 for different water vapor concentrations 
indicate that the radiative interaction increases slowly with an increase in the amount of the gas. 
The results for 50% H 2 O are presented in Fig. 5 for two different pressures (Poo=l ar >d 3 atm) 
and x-locations (x=5 and 10 cm). It is noted that the increase in pressure has dramatic effects 
on the radiative interaction. 

For a mixture of 50% H 2 O in air, the conduction and radiation heat transfer results are 
compared in Fig. 6 for Poo = 3 atm and for two different x-locations (x=5 and 10 cm). The results 
demonstrate that the conduction heat transfer is restricted to the region near the boundaries and 
does not change significantly from one x-location to another. The radiative interaction, however, 
is seen to be important everywhere in the channel, and this can have a significant influence on 


the entire flowfield. 
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Comparative results for 100% homogeneous species of H 2 0, NO and OH are illustrated in 
Fig. 7 at the exit plane (x=10 cm). As would be expected, the radiative contribution of H 2 0 (with 
five bands) is significantly higher than NO and OH. Only the fundamental bands of NO and OH 
are considered in this study, and it is noted that NO is a better radiating gas in comparison to OH. 

For a mixture of 25% H 2 0 in air, radiative flux for two different plate spacings (L=3 and 
6 cm) are illustrated in Fig. 8 for two x-locations (x=5 and 10 cm). The rate of radiative 
transfer is a strong function of the amount of the participating species and the pressure path 
length, PL. consequently, the results for the larger plate spacing indicate significantly higher 

radiative interactions. 

The effect of increased Mach number on the radiative transfer is illustrated in Fig. 9 for 
pure H 2 0 and for a mixture of 50% H 2 0 in air. The results shown are for the exit plane (x=10 
cm). At higher Mach number, the boundary layer is relatively thinner and the temperature in the 
boundary layer is significantly higher. This, in turn, results in higher rate of radiative transfer. 

For a mixture of 50% H 2 0 in air, comparative results for the parallel plate channel and the 
circular tube are presented in Fig. 10 for two x-locations (x=5 and 10 cm). The results for the 
circular tube in general, exhibit the same trend as for the parallel plate geometry. Since the 
circular geometry provides additional degrees of freedom for radiative interactions [17], the rate 
of radiative transfer is higher for the tube. 

The influence of radiative interactions in chemically reacting supersonic internal flows 
was investigated by considering the physical model shown in Fig. lc. The specific problem 
considered is the supersonic flow of premixed hydrogen and air in a channel with a compression 
corner on the lower boundary. The physical dimensions considered for obtaining results are 
L=2 cm, x 1= 1 cm, x 2 =2 cm, L*=xi+x 2 =3 cm, and a=10°. The inlet conditions, which are 
representative of the scramjet operating condition, are P=1 atm, T=900 K and M=4.0. The 
flow is ignited by the shock from the compression comer. The flowfield for this problem has 
been investigated by several researchers [21-28] where different chemistry models have been 
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used. Influence of radiative interactions was investigated in [25-27] by considering a simple two- 
step chemistry model. Recently [28], a comparative study of the flowfield was conducted by 
employing three chemistry models (see Table 1). It was found that significant amount of radiation 
participating species are produced by the 35-step chemistry model. The results presented in Figs. 
3-10 provided essential information on radiative behavior of these species. Complete discussions 
on use of the three chemistry models with and without radiation are provided in [28]. Selected 
results are presented here to demonstrate the influence of radiative interactions. 

The computed results for the 35-step chemistry model are presented in Figs. 1 1-14 with and 
without radiative interactions. The results were obtained by using a 31x31 grid; this was found 
to be an appropriate grid for the model. The variations in temperature, pressure, and species 
concentrations along the x-coordinate are shown for a y-location of 0.02 cm from the lower wall. 
The shock occurs at about x/L*=l/3, and it is noted that the 35-step chemistry model predicts 
the ignition time accurately. The temperature is seen to increase uniformly along the channel 
(Fig. 11) and there is a significant increase in pressure after the shock (Fig. 12). Results of 
Figs. 13 and 14 show that significant amounts of radiation participating species H 2 O and OH are 
produced after the shock. The effect of radiative interaction is to lower the amount of species 
production due to radiative transfer in the x-direction. 

The results of radiative transfer by the three chemistry models are compared in Figs. 15 and 
16 at y=0.02 cm. The results for the normal radiative flux presented in Fig. 15 demonstrate that 
radiative interactions increase rapidly after the shock. The three models are seen to predict the 
same general trend. The results of streamwise radiative flux illustrated in Fig. 16 show that the 
net q Rx decreases towards the end of the channel. This is due to cancellation of fluxes in the 
positive and negative x-directions. It is noted that the net radiative transfer is in the negative 
x-direction. The 18-step and 35-step models are seen to predict significantly higher qRx than 
the 2-step model. This is because radiative heat transfer is a strong function of temperature, 
pressure and species concentration which are higher (in the positive x-direction) for the 18-step 
and 35-step models than the 2-step model. 
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Conclusions 

Two-dimensional compressible Navier-Stokes equations have been used to investigate the 
influence of radiative energy transfer on the entrance region flow under supersonic flow condi- 
tions. Computational procedures have been developed to incorporate gray as well as nongray 
formulations for radiative flux in the general governing equations. Specific results have been 
obtained for different amounts of H 2 O, OH and NO in combination with air. Results demon- 
strate that the radiative interaction increases with an increase in pressure, temperature and the 
amount of participating species. This can have a significant influence on the overall energy 
transfer in the system. Most energy, however, is transferred by convection in the flow direction. 
The radiative interactions in reacting flows have been investigated by considering the supersonic 
flow of premixed hydrogen and air in a channel with a compression comer at the lower bound- 
ary. Depending upon the chemistry model employed, the radiative interaction is seen to change 
significantly the temperature, pressure and species concentration in the flow direction. 
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Tfcbte 1. Hydrogen-Air Combustion Mechanism (40) 


REACTION 

A(moles) 

N(cm J ) 

E(cak>r ics/gn i - mole ) 

•• following reactions consulate the 18-step model •• 


(1) Oi + II) «-* OH + OH 

1.70 * 10'* 

0 

48150 

(2) Oj + II «-• OH ♦ O 

1.42 x I0' 4 

0 

16400 

(3) H) + OH +. HjO + H 

3.16 x HI 7 

1.8 

3030 

(4) Hi + O «-► OH + H 

2.07 x HI 14 

0 

13750 

(5) OH ♦ OH 11)0 + O 

5.50 x 10'* 

0 

7000 

(6) HtOiltMH 11)0 + M 

2 21 x 10“ 

—20 

0 

(7) II + II ♦ M « • II) + M 

6.53 x I0' 7 

— 1.0 

0 

(H) II + Oi + M « » IK >2 ♦ M 

3.20 x 10" 

10 

0 

(9) OH + HO) »♦ O) ♦ 11)0 

5.IX) x I0 11 

0 

KMX) 

(10) II + IIO) «-» II) + O) 

2.53 x 10'* 

0 

7fXJ 

(11) 11 + IIO) «-* OH ♦ OH 

1.99 x I0 14 

0 

1800 

(12) O + HO) *-» O) -*• OH 

5.00 x III" 

0 

KXXJ 

(13) HQj(HO) *-*0) + 11)0) 

1.99 x 10 12 

0 

0 

(14) II) -r IIO) ++ H + H)0) 

3.01 x 10" 

0 

18700 

(13) OII + 11)0) ++ H)0 + HO) 

1.02 x 10" 

0 

1900 

(16) II + 11)0) ++ H)0 + OH 

5.00 x I0 14 

0 

10000 

(17) O + 11)0) ++ OH + HO) 

1.99 x 10° 

0 

5900 

(18) 11)0) + M ++ OH + OII + M 

1.21 x I0 17 

0 

45 SIX) 

♦♦ remaining reactions complete Uie 35 -step model •• 


(19) O) + M +♦ O ♦ O + M 

2.75 x 10" 

—10 

118700 

(20) N] + M«-»NtN + M 

3.70 x 10*' 

—1.6 

225000 

(21) N + O) O + NO 

6.40 x 10* 

1.0 

6300 

(22) N ♦ NO ++ O ♦ N) 

1 60 x 10* J 

0 

0 

(23) N + Oil «-* II ♦ NO 

6.30 x 10“ 

05 

0 

(24) 11 + NO + M +-> UNO + M 

5.40 x I0 15 

0 

—600 

(23) II + UNO «-» II) + NO 

4 80 x III" 

0 

0 

(26) O + UNO «-» OH + NO 

500 x 10" 

0.5 

0 

(27) OH + UNO +♦ H)0 + NO 

3 60 x 10" 

0 

0 

(28) IIO) + UNO +* 11)0) + NO 

2.00 x 10" 

0 

0 

(29) HOi + NO ++ OH + NOj 

3.43 x 10" 

0 

— 260 

(30) II + NO] OH + NO 

3.50 x 10 14 

0 

1500 

(31) 0 + NO) n O) + NO 

1.00 x 10" 

0 

6<K) 

(32) NO) + M «-» O + NO + M 

1.16 x 10" 

0 

66000 

(33) M + OH + NO ++ UNO) + M 

5.60 x 10" 

0 

— 1700 

(34) M 4- Oil ♦ NOj ♦-# HNOj + M 3.00 * 10 ,s 

0 

—3800 

(35) OH + UNO) +♦ 11)0 + NO) 

1.60 x I0 i2 

• 

0 

0 

*♦ following reactions constitute the global 2-step model 124. 27 , 39] ♦* 

(1") H] ♦ O] «-» 2 OH 

11.4 x HI 47 

— 10.0 

4865 

(2") 2 OH + H) *-» 2 H)0 

2.50 x I0 64 

— 13.0 

42500 1 





Fig. la Radiating gas flow between parallai boundaries. 



Fig. lb Radiating gas flow within a circular tube. 



Fig. lc Radiating gas flow in a channel with a compression 
expansion ramp. 
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Fig. 2 Divergence of rsdive flux along the channel for gray 
and nongray models, M. » 4.3. 
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Fig. 4 Radiativi flu* vs. y at the channel exit, M-, ■ 3.0. 


Fig. 6 Radiative and conductive fluxee vs. y for P v3 atm, x * 5 
and 10 cm, 60% H 2 0, ML * 3.0. 
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Fig 6 Radiative flux vt. y for P ■ 1 and 3 atm, i » 6 and 10 cm, 
60* H*0, M_ . 3.0. 
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Fig. 7 Radiative flu* vs. y for HjO, NO, end HO, P - 1 etm. 
M_ - 3.0. 
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Fig. 8 Radiative flux vs. y for two different plate epeerngs (L ■ 3 
and 6 cm) at x » 5 and 10 cm, M_ * 3.0. 


Fig. 10 Radiative flux vs y for channel and tube (L » D « 3 cm) 
at x • 6 and 10 cm, 50% HiO, M. * 3.0. 
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Fig. 9 Radiative flux re. y for 50% and 100% HjO, x « 10 cm, 
• 3.0 and 4.6. 


Fig. 11 Temperature variation with x for reacting, and reacting 
and radiating flows. 
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Pressure ratio , P / P 4 



Fig. 12 P re mu re variation with x for reacting, and reacting and 
radiating flows. 



Fig. 13 Variation of H 2 O maw fraction with x for reacting, and 
reacting and radiating flows. 



Fig. 14 Variation of OH mass fraction with x for reacting, and 
reacting and radiating flows. 
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